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We study the algebraic structure of the Dirac fermion state in the organic conductor, a- 
(BEDT-TTF) 2 i3. We find a pair of generators for the Hamiltonian of a-(BEDT-TTF) 2 l3 that 
describes the chirality of Dirac fermions. The phase parameters associated with those genera- 
tors have an intimate and simple relationship with the positions of the Dirac points in the Bril- 
louin zone. By making use of the form of the generators, a reduced form of the Hamiltonian 
is constructed that shares some characteristic features of the Dirac fermions with ar-(BEDT- 
TTF) 2 i3. For the reduced Hamiltonian we present the analytic expression for the Dirac point 
position that sits at arbitrary point in the Brillouin zone determined by the hopping parame- 
ters. 
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1. Introduction 

Massless Dirac fermion spectrums in two-spatial dimension are now realized in vari- 
ous condensed matter systems. A well-known system is graphene, a single atomic layer of 
graphite, which is a two-dimensional hexagonal lattice of carbon atoms. 1 The existence of the 
Dirac fermion spectrum was clearly demonstrated by the observation of the half-integer quan- 
tum Hall effect. 2 3 Another example is a group of three-dimensional topological insulators 4-6 
where the surface state exhibits a Dirac fermion spectrum. A Dirac fermion spectrum appears 
even in superconductors. In the antiferromagnetic state of iron pnictide superconductors, a 
Dirac fermion spectrum 7,8 was observed in angle-resolved photoemission spectroscopy. 9 An 
organic compound, a-(BEDT-TTF) 2 I 3 [BEDT-TTF=bis(ethylenedithio)tetrathiafulvalene], 
which has a layered structure of BEDT-TTF molecules, demonstrates that a Dirac fermion 
spectrum appears by suppressing strong electron correlation effects. 10 Under ambient pres- 
sure, this compound undergoes a metal-insulator transition due to a charge ordering induced 
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by a strong electron correlation. 1114 This charge ordering is suppressed by applying pres- 
sure, and transport measurements support the presence of a Dirac fermion spectrum 1518 for 
pressures higher than 1 .5GPa. The presence of Dirac fermions is clearly demonstrated in the 
interlayer magnetoresistance measurement 19 where the zero energy Landau level of Dirac 
fermions leads to negative magnetoresistance. 20 

A natural question is how a Dirac fermion spectrum appears. In case of graphene, the 
mechanism is trivial. The Dirac fermion spectrum is realized as a result of symmetry of the 
honeycomb lattice. In case of topological insulators, essentially the spin-orbit interaction lead 
to a linear Dirac fermion energy spectrum. 4 6 ' 2123 In case of iron-pnictide superconductors, a 
characteristic configuration of Fe and As leads to degeneracy of d-orbitals at the Y point, and 
supports the appearance of Dirac fermion spectrum. 7 Contrary to these systems, the mecha- 
nism of stabilizing the Dirac fermion spectrum in a-(BEDT-TTF) 2 I 3 is elusive. 24 ' 25 

In Dirac fermion systems the conduction and the valence bands touch at two inequiva- 
lent points in the Brillouin zone (BZ). Usually reflecting symmetry of the system, the Dirac 
points appear at high symmetric points or along symmetric lines in the BZ. However, the 
Dirac points in a-(BEDT-TTF) 2 l3 are neither located at symmetric points of the BZ nor along 
symmetric lines in the BZ. So the presence of the Dirac fermion spectrum of this system 
is called accidental degeneracy. 26 Furthermore, the positions of the Dirac points depend on 
the applied pressure. 10 ' 27 The complicated four band structure of aKBEDT-TTF)^ prevents 
revealing the essential algebraic structure for the stabilization of the Dirac fermion spectrum. 

In this paper, we study the underlying algebraic structure for the Dirac fermion spectrum 
in a-(BEDT-TTF) 2 l3. In spite of the complexity of the system, i.e., the presence of many 
hopping parameters and multi-orbital character, there is a relatively simple algebraic struc- 
ture that supports the Dirac fermion spectrum. We show that two matrices Y x and Y y describe 
chirality of Dirac fermions, and, interestingly, phase parameters contained in these matrices 
and the positions of the Dirac points are intimately connected in a simple functional rela- 
tionship. Using the result, we construct a multi-orbital Hamiltonian with a Dirac fermion 
spectrum whose Dirac points move in the BZ by changing the hopping parameters. 

The paper is organized as follows. In § 2 we consider graphene as a simple example to 
illustrate the algebraic structure of a Dirac fermion system. In § 3 we introduce the model for 
a-(BEDT-TTF) 2 I 3 , and we show a simple relationship between the matrices Y x and Y y and the 
positions of the Dirac points. In § 4 we construct a multi-orbital Hamiltonian that has a Dirac 
fermion spectrum whose Dirac points move by changing the hopping parameters. We also 
present the analytic expression for the Dirac point position in the BZ. In § 5 we summarize 
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the result. 



2. Algebraic structure of graphene 

In order to illustrate the algebraic structure of Dirac fermions we consider graphene as 
a simple example. The purpose is to extract the algebraic structure that is associated with 
chirality, which is one of the most characteristic properties of Dirac fermions. 

Usually chirality is defined by using eigenstates of Dirac fermions. 28 In case of two 
component Dirac fermion wave functions, which is the most simplest Dirac fermion wave 
function, one can define a chirality vector field by n ka = (^k\c a \i//k) with if/ k being a two- 
component Dirac fermion wave function at k in the Brillouin zone and <j a being the Pauli 
matrices. The chirality vector has either a vortex or anti- vortex configuration around Dirac 
points. It is possible to construct this chirality vector field directory from the Hamiltonian by 
introducing two matrices or generators. We illustrate how this construction is carried out in 
graphene. 

The Hamiltonian for graphene is given by 28 



c Ak L Bk 



Hi 



( \ 

CAk 
CBk ) 



(1) 



where c k (c ak ) creates (annihilates) an electron with the wave vector k on sublattice a = A, B, 



and 



H k = 





* 



Kk 





with 



K k = -t 



^ + e {-\h + fk y ) + e -{\k x+ 4k y ) 



(2) 



(3) 



We set the lattice constant unity. Denoting a Dirac point in the BZ as k D , we find K ko = 0. 
Expanding K k around k = k D , we find that 



dk x 



k=k f 



= C x exp {-^Dx^j i, -ttt =C y exp f-^o*) , 

y k=k D 



where C x and C y are real constants. Note that there is the factor exp(-^fc Dv ) that depends 
on k/) in the right hand side of these equations. We encounter similar factors for the case of 
a-(BEDT-TTF) 2 l3 later. Introducing the following matrix, 

^ dH k 



dk a 



k=k c 



and representing these matrices as Y x = a x (cr x cos 9 X + cr y sin 6 X ) and r v = a y (cr x cos 9 y + 
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(T y sin#y), we find \9 X - 9 y \ = n/2. This phase difference n/2 comes from the phase difference 
between dK k /dk x \ k=k and dK k /dkJ . If we define an inner-product between two matrices 

u - Ik— 

as (A, B) = tr(A5) = JjjjAjjBjj, we find that (T x , T y ) = a x a y cos(6 x - 6 y ) = 0. because of the 
phase difference n/2. 

Now we construct the chirality vector directly from the Hamiltonian using the matrices 
r a . The 2x2 matrix H k has the following form around the Dirac point: 

H kD+sk = T x 5k x + T y 8k y + ■■■ . (4) 

Suppose one moves around the Dirac point. The path is represented by Sk x = Sk cos(p, 
Sk y = 6k sin (p, with (p changing from to 2n. We define a vector field (n kx , n^) as 

n kx = ix{T x H k ), n ky = tr(r y H k ), (5) 

It is easy to see that n kx and n^, vanish at Dirac points. As shown in Fig. 1, we find that this 
vector field rotates as one goes around the Dirac point. The vector field (n kx , n ky ) forms vortex 
(anti- vortex) around the two inequivalent points. These points correspond to the Dirac points 
in graphene. Therefore, the vector field (n kx , n ky ) describes the chirality of the Dirac fermion. 
Note that by symmetry the vector field (n kx , n^) can be rotated in the k x - k y plane, or one can 
define 

n' kx = tr(Y' x H k ), n' ky = tr(V y H k ), (6) 

with T' Y = r Y cos (p - r v sin (p and Y' Y = Y x sin (p + T y cos (p. In Fig. 1 we set T x = cr x and r y = cr y 
by making use of this symmetry. 

In case of graphene, the analysis is trivial because we have only four matrices, the Pauli 
matrices cr x , cr y , cr z , and the unit matrix cr , to construct 2x2 Hamiltonians. Because of the 
absence of <x and cr z components in the graphene Hamiltonian, the space of matrices T a is 
uniquely determined. 

3. Algebraic structure of (r-(BEDT-TTF) 2 I 3 

Contrary to graphene, it is unclear how to find T a matrices for the case of ff-(BEDT- 
TTF) 2 l3 because the Hamiltonian is 4 x 4. (See, eq. (8) below.) There are 15 generators 
and one unit matrix for constructing 4x4 Hermite matrix. A suitably chosen set of those 
generators satisfy commutation relations similar to that for cr x , cr y , and o%. However, naive 
choices fail to describe chirality of Dirac fermions in a-(BEDT-TTF)2l3. 

In order to find Y a , we study the Hamiltonian of a-(BEDT-TTF) 2 l3. We consider a single 
plane of BEDT-TTF molecules. The unit cell consists of four BEDT-TTF molecules, which 
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Fig. 1. The vector field (rik X , «jt v ) for graphene in k space. The vector field is normalized to clarify the structure 
around the Dirac points. (We take the lattice constant as the unit of length.) The vector field forms vortex (anti- 
vortex) around K:(2^/3)(l, 1 / V3) (K':(2;r/3)(1, -1/ V3)). The points K and K' are the Dirac points. (A similar 
result was presented in ref. 29.) 



are labeled by 1,2, 3, and 4. The Hamiltonian is given by 10 30 

( H = J^clH k c k , 



(7) 



where the spinor representation of the creation operator is c 



C lk °2k C 3k C 4k 



and 





t c i + t cZ e 



iky 



t c i + t c2 e 




-iky 



t. A e 



t pl -t p4 e- ik * t p4 e~ ik y -t pl e< k ^ 



. t P 2 - t p3 e ik * t p3 - t p2 e ,K * t c3 + t c4 e" 

Here hopping parameters t n (rj = cl, ...,c4, pi, ...,p4) are assumed to be pressure dependent. 
There are two types of pressure: uniaxial pressure and hydrostatic pressure. For the purpose 
of investigating algebraic structure of eq. (8) in wide hopping parameter ranges, we consider 
a uniaxial pressure applied along the a-axis, which is the stacking direction of BEDT-TTF 
molecules in the plane. For the pressure dependence of t v = tr,(P a ), we take extrapolation 
formula given by 30 



-a. 



t p i — tp\e 
ik y - t^eK^y) 


Jky 



t P 2 — t p ?,e 
t P 3 - t p 2e 
t C 3 + t C Ae 




ik x 



iky 



-iky 



(8) 



t v (Pa) = t„(0)(l + K n P a ). 



(9) 
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The parameters ^(0) and K n are determined from the data at P a = given in ref . 3 1 and the 
data at P a = 2kbar given in ref. 32. The explicit values are 30 t pl (0) = 0.140, t p2 (0) = 0.123, 
t p3 (0) = -0.025, t p4 (0) = -0.062, f cl (0) = 0.048, t c2 (0) = -0.020, t c3 (0) = -0.028, t c4 (0) = 
-0.028, K pl = 0.011, K p2 = 0, K p3 = 0, K p4 = 0.032, K cl = 0.167, K c2 = -0.025, and 
K c3 = K c4 = 0.089. 

Our purpose is to find two matrices T x and Y y that play a role similar to cr x and cr y in 
graphene. Contrary to the Hamiltonian of graphene (1), there are other orbitals. Therefore, 
the situation is much more complicated. In fact, one can find Y x and Y y from numerical cal- 
culations at a given pressure P a . However, it is hard to figure out the underlying algebraic 
structure. Nevertheless we show that there is a simple algebraic structure that governs the 
Dirac fermion spectrum of Q'-(BEDT-TTF) 2 l3. We assume that T x and T y have the following 
form: 























e ie 


e -i4>x 














e -i9 x 
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g-Upy 




















e w v 








e -i9y 






(10) 



These forms are inferred from the forms of dHk/dk x and dH^/dky. Some components are 
omitted so that tr^Fy) = 0. Furthermore, we take Y x and Y y so that n^ x = tx(Y x Hk) depends 
on tp\,t p2 , t P 3, t p4 and n ky = ir(Y y H k ) depends on t c i,t c2 , t c3 , t c4 . 

If we represent a Dirac point by (k Dx ,k Dy ), the phases in eq. (10) should depend on k Dx 
and koy In order to make clear this dependence, we shift 6 X and <p y as 6 X —* 6 X + k Dx /2 and 
<p y —* <p y - k Dy /2, respectively, and fix cf) x and 9 y as <p x = n/2 + k Dx /2 and 9 y = n/2 - k Dy /2, 
respectively. We find that the vector field {n kx , nky) has a vortex configuration in k space as 
shown in Fig. 2 on the plane of k x and k y . The vortex moves in k space by changing the phase 
parameters <f> y and 9 X . The chirality of Dirac fermions in a-(BEDT-TTF)2l3 is described by 
("fa, nky) by suitably choosing the phase parameters. 

The phase parameters are determined by solving (nk DX , nk D} ) = 0: 



. . k Dx I k Dx \ I k Dx \ . k Dx 

-t p i sin — - t p2 cos 1 6 X — — I + t p3 cos \9 X + — I - t p4 sin — = 0, 



(11) 



k D y\ (. k Dy \ k Dy k Dy 

t cX cos | <f> y - — I + t c2 cos I <p y + — I + t c3 sin — - t c4 sm — = 0. 



(12) 



The pressure dependence of the phase parameters is shown in Fig. 3. The Dirac points k D 
are computed by diagonalizing the Hamiltonian numerically. Figure 4 shows the pressure 
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Fig. 2. The vector field (n**, n^) for a-(BEDT-TTF)2l3 in k space at Pa = 5kbar. The vector field is normal- 
ized to clarify the structure around the Dirac point as in Fig. 1. The phase parameters are 9 X = 161.85 degree, 
</> y = 50.05 degree. 



dependence of the positions of the Dirac points k D . Note that the pressure dependence of 9 X 

is similar to that of k Dx and the pressure dependence of cp y is similar to that of k Dy . Figure 5 

shows k Dx dependence of 6 X and k Dy dependence of <f> y . We find that the phase parameters are 

approximately linearly dependent on the Dirac point coordinates, that is, 

2/ 3 \ 
Ox = ^ \ kox + 4 n j ' 

cp y = | {k Dy + . 

This result demonstrates that underlying symmetry structure is rather simple despite the fact 
that the Hamiltonian has a complex form. The reason why we encounter the factor 2/3 has 
not been clarified yet. The origin of this simple factor is left for a future study. 

4. Multiorbital Hamiltonian with Dirac fermion spectrum 

In the previous section, we have found matrices Y x and Y y that describe chirality of 
the Dirac fermions in Qr-(BEDT-TTF) 2 i3. In this section, we construct a multi-orbital Dirac 
fermion Hamiltonian where chirality is exactly described by T x and T y . We show that the 
Hamiltonian has Dirac points whose positions are hopping parameter dependent and not lo- 
cated at symmetric points in the BZ. So characteristic features of Dirac fermions in a-(BEDT- 
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Fig. 3. Pressure dependence of the phase parameters 6 X and tp y . 
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Fig. 4. Pressure dependence of the Dirac point in the third quadrant of BZ. 



TTF) 2 l3 are reproduced by this Hamiltonian. 

From the form of eq. (10) and the fact that the chirality vector field is obtained by taking 
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Fig. 5. The phase parameters 9 X and tf> y as functions of Hdx and ko y , respectively. As shown by a dashed line, 
6 X and <j> y approximately linearly depend on k Dx and k Dy , respectively. 



trace of T x H k and T y H k , we may consider the following Hamiltonian: 






p + qe lky 


r + se lkx 





p + qe lky 








f + ue ikx 


r + se~ lkx 








v + we~ iky 





t + ue~ ikx 


v + we lky 






The hopping parameters, p, q, r, s, t, u, v, and w, which are real numbers, are defined in Fig. 6 
in real space. The Hamiltonian (13) is exactly diagonalized. The eigenvalues are 

E { r ) = ±^A k ±^^ k (14) 

where 

A k = + qe~ iky \ 2 + \r + se ikx f + \t + ue ikx f + \v + we~ iky f\ , (15) 

B k = \(p + qe' iky )(v + we ikv )-(r + se ikx )(t + ue~ ikx )\. (16) 

The upper two bands E ( k + ' +) and EY'~' touches at a point when A k = B k is satisfied at 
a single k point. Since there are many possibilities to satisfy this condition in general, we 
limit ourselves to consider a simple case where the equations v + we~' ky = e'^ (p + qe~ lk A and 
t + ue lkx = -e'^ (r + se lkx \ are satisfied at some point (k x ,k y ) with <p being a constant. These 
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Fig. 6. The hopping parameters of the reduced Hamiltonian eq. (13). The unit cell contains four sites as in the 
case of a-(BEDT-TTF) 2 I 3 . 



equations are easily solved analytically with respect to k x and lc y , and we find 

t + re i4 



k r = Im 



In 



k Y = -Im 



In 



u + se l( ? 
pf M> _ v 



K w - qe" 

with constraints \pe"^ - v\ = \w - qe x<t> \ and \t + re'^l = \u + se l % The constraints are solved with 
respect to w and s: 

w = q cos <p± yjp 2 + v 2 - 2pv cos (p - g 2 sin 2 0, (17) 



s = u cos <p + yjr 2 + t 2 + 2rt cos <p - w 2 sin 2 </>. 



(18) 



We show an example in Fig. 7 which always exhibits a zero gap state, i.e., the energy of the 
Dirac point does not intersect either E^' +) or Erf' _) in the BZ except at the Dirac points. We 
clearly see that the energy dispersion is linear around {k x ,k y ) = (0.3587, 1.773)^, which is 
the Dirac point. Although we do not know the analytic solution for a-(BEDT-TTF) 2 i3 but 
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Fig. 7. The energy dispersions of the Hamiltonian (13) along k x -axis with k y = 1.773tt (left panel) and £ v -axis 
with k x = 0.35877T (right panel). The upper two bands, i.e., and £1 are shown. The Dirac point is at 

(k x , k y ) = (0.3587, 1 .773)7r. The horizontal solid line is the energy of the Dirac point. The parameter values are 
p = 0.08808, q = -0.0175, r = 0.1477, s = -0.0389, t = -0.025, u = -0.123, v = -0.04046, w = 0.1015, 
and <p = 70 degree. We set these values using the hopping parameter values of ff-fBEDT-TTF^fi at P a = 5kbar 
except for s and w. The parameters s and w are determined from eqs. (17) and (18). 



we believe that the Hamiltonian (13) captures essential symmetric properties of the Dirac 
fermions of a-(BEDT-TTF) 2 i3. 

The fact that the Dirac point moves in the BZ by changing the hopping parameters is 
demonstrated in Fig. 8. Therefore, the Hamiltonian eq. (13) shares the feature of the moving 
Dirac points in the BZ with a-(BEDT-TTF) 2 l3. We show one Dirac point k D = (k Dx ,k Dy ) in 
Fig. 8. The other Dirac point is at (-k Dx , -k Dy ) by symmetry. 

It is interesting to point out a possible relationship between the model and the ^-flux state 
that was originally proposed for the two-dimensional antiferromagnetic Heisenberg model. 33 
We shall discuss that a-(BEDT-TTF) 2 l3 is characterized by a 7r-flux state in a multi-orbital 
system in a future publication. 

5. Summary 

To summarize, we have investigated underlying algebraic structure of the Dirac fermion 
spectrum of a-(BEDT-TTF) 2 l3. We have found that chirality of Dirac fermions in a-(BEDT- 
TTF) 2 I 3 is well described by the matrices Y x and Y y defined by eq. (10). The phase factors in 
T x and T y have intimate relationships with the position of the Dirac point in the BZ, and we 
have found a simple relationship between them. A reduced form of the Hamiltonian has been 
constructed from Y x and T y . The Dirac point of the reduced Hamiltonian moves in the BZ 
by changing hopping parameters. Although the analytic expression for the Dirac points of a- 
(BEDT-TTF) 2 I 3 is still unknown, our reduced Hamiltonian has the analytic expression for the 
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Dirac points. Therefore, the reduced Hamiltonian is a useful model to unveil the mechanism 
of stabilizing the Dirac fermion spectrum in ar-(BEDT-TTF) 2 i3. 
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